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Abstract. Consider a N x n non-centered matrix En with a separable variance profile: 

■ An . 



Matrices Dn and Dn arc non-negative deterministic diagonal, while matrix An is de- 
terministic, and Xn is a random matrix with complex independent and identically dis- 
tributed random variables, each with mean zero and variance one. Denote by Qn{z) the 
resolvent associated to S 

Qn{z) = (s„s; - 2/jv)-l . 

Given two sequences of deterministic vectors iun) and with bounded Euclidean 

norms, we study the limiting behavior of the random bilinear form: 

U*„Q„{z)Vn , 6C-M+ , 

as the dimensions of matrix E„ go to infinity at the same pace. Such quantities arise in 
the study of functionals of S„E* which do not only depend on the eigenvalues of , 
and are pivotal in the study of problems related to non-centered Gram matrices such as 
central limit theorems, individual entries of the resolvent, and eigenvalue separation. 
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1. Introduction 

The model. Consider a N x n random matrix I]„ ~ (^j") given by: 

1 ^ j_ 

= 7= I" = Yn -\- A„ , (1-1) 

Vn 

where D„ and Dn are respectively N x N and n x n non-negative deterministic diagonal 
matrices. The entries of matrices (X„), (X^" ; n) are complex, independent and identically 
distributed (i.i.d.) with mean and variance 1, and An — (afj) is a deterministic N x n 
matrix whose spectral norm is bounded in n. 

The purpose of this article is to study bilinear forms based on the resolvent Qn{z) of 
matrix I]„I]* , where S* stands for the hcrmitian adjoint of 



Qn{z) = (S„i;* - zIn) ^ 
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as the dimensions and n grow to infinity at the same pace, that is: 

. . N . N 

< hminf — < hmsup — < oo , (1-2) 
n n 

a condition that will be referred to as A^, n — > cx) in the sequel. 

A lot of attention has been devoted to the study of quadratic forms y* Ay, where y = 
• • ■ XnY' , the XiS being i.i.d., and A is a matrix independent from y. It is well- 
known, at least since Marcenko and Pastur's seminal paper [151 Lemma 1] (see also [H 
Lemma 2.7]) that under fairly general conditions, y* Ay 7i~^Tr A. 

Such a result is of constant use in the study of centered random matrices, as it allows to 
describe the behavior of the Stieltjes transform associated to the spectral measure (empirical 
distribution of the eigenvalues) of the matrix under investigation, see for instance [23], [24], 
[14l[T5], etc. Indeed, the Stieltjes transform of the spectral measure writes: 

1 1 ^ 

Uz) = ^TrQ„(z) = -Y}Qn{z)Uz) , 

1=1 

where the [Qn{z)]iiS denote the diagonal elements of the resolvent. Denote by fji the ith 
row of Yin and by matrix I]„ when row fji has been removed, then the matrix inversion 
lemma yields the following expression: 

In the case where S]„ — n^^^^Xn, the quadratic form that appears in the previous expression 
can be handled by the aforementioned results. However, if I]„ is non-centered and given by 
(jl.ip . then the quadratic form writes: 

ViQi{z)fl* = mQi{z)yl + a^Qi{z)y* -f mQ,{z)a* + aiQ,{z)a* , 

where Qi{z) = (S* jl]„,i — zl)~^, and iji and di are the ith rows of matrices Yn and An- 
The first term can be handled as in the centered case, the second and third terms go to 
zero; however, the fourth term involves a quadratic form aiQi{z)d* based on deterministic 
vectors. 

It is of interest to notice that, due to some fortunate cancellation, the particular study 
of bilinear forms of the type UnQn{z)vn or their analogues of the type u„Q„(z)i;* can be 
circumvented to establish first order results for non-centered random matrices (see for in- 
stance [8], [15]). However, such a study has to be addressed for finer questions such as: 
Asymptotic behavior of individual entries of the resolvent (see for instance [TUl Eq. (2.16)] 
where such properties are established in the centered Wigner case to describe fine properties 
of the spectrum) , Central Limit Theorems [TT] [13] , behavior of the extreme eigenvalues of 
, behavior of the eigenvalues and eigenvectors associated with finite rank perturbations 
of SnS* [6], behavior of eigenvectors or projectors on eigenspaces of Q{z) (see for instance 
[3] in the context of sample covariance (centered) model), etc. 

In a more applied setting, functionals based on individual entries of the resolvent [T] 
naturally arise in the field of wireless communication (see for instance Section IO]) . Moreover, 
the asymptotic study of the quadratic forms m*(5„(z)w„ is important in statistical inference 
problems. In the non-correlated case (where Z)„ = /jv and D„ = /„), it is proved in |25| how 
such quadratic forms yield consistent estimates of projectors on the subspace orthogonal to 
the column space of An in the Gaussian case (see also Section [^?^ . Such projectors form the 
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basis of MUSIC algorithm, very popular in the field of antenna array processing. A similar 
approach has been developed in jl9| . |20j for sample covariancc matrix models. 

It is the purpose of this article to provide a quantitative description of the limiting behavior 
of the bilinear form u*Q„(z)?;„, where u„ and w„ are deterministic, as the dimensions of I]„ 
go to infinity as indicated in (|1.2p . 

Assumptions, fundamental equations, deterministic equivalents. Formal assump- 
tions for the model arc stated below, where || • || either denotes the Euclidean norm of a 
vector or the spectral norm of a matrix. 

Assumption A-1. The random variables {X^- ; l<ii<N,l<i<n,n>l) are complex, 
independent and identically distributed. They satisfy EX^" = and EjXJjP = 1. 

Assumption A-2. The family of deterministic N x n matrices > 1) is bounded for 

the spectral norm as N,n oo: 

dmax = sup ||A„|| < OO . 
n>l 

Notice that this assumption implies in particular that the Euclidean norm of any row or 
column of \\An\\ is uniformly bounded in N,n. 

Assumption A-3. The families of real deterministic N x N and n x n matrices (Dn) and 
(Dn) are diagonal with non-negative diagonal elements, and are bounded for the spectral 
norm as N,n —i- cg: 



Moreover, 



dmax = sup II Aill < OO and dmax = SUp || Ai|| < OO . 
n>l n>l 



dmin = inf —TrDn > and dmin = inf — Tr_D„ > 

N N n n 



We collect here results from [15]. 
The following system of equations: 



z e 



5{z) = iTri?„ (^-z(/Ar + ^(z)i?„)/Ar + A„(/„ + <5(z)^„) 

(1.3) 

admits a unique solution [5, 5) in the class of Stieltjes transforms of nonnegative measure^ 
with support in R"*". Matrices T„(z) and T„(z) defined by 



r„(z) = l^-z{lN + 5{z)an) + A,,(ln + 5{z)D^^ ^ A;)j 

fn{z) = l^-z{In+5{z)bn)+ A*^(In +~6{z)Dn) ' 



(1.4) 



^In fact, S and & axe the Stieltjes transforms of measures with respective total mass n ^TrZ)„ and 
n-iTrD„. 
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are approximations of the resolvent Qn(z) and the co-resolvent Qn{z) — (S* S„ — ^/^v) ^ in 
the sense that (JtdL^ stands for the almost sure convergence): 

i-Tr(Q„(z)-r„(^))-^^^^0, 

which readily gives a deterministic approximation of the Stieltjes transform N^^Ti Q„ {z) of 
the spectral measure of SnSJj in terms of T„ (and similarly for Qn and r„). Matrices T„ 
and Tn will play a fundamental role in the sequel. 

Nice constants and nice polynomials. By nice constants, we mean positive constants 
which depend upon the limiting quantities dmin, dmin, <imax, dmax, Omax, liminf and 
limsup ^ but are independent from n and N. Nice polynomials are polynomials with fixed 
degree (which is a nice constant) and with non-negative coefficients, each of them being a 
nice constant. Further dependencies are indicated if needed. 

Statement of the main result. Let dz be the distance between the point z G C and the 
real nonnegative axis M^: 

Sz ^ dist(z,]R+) . (1.5) 
Here is the main result of the paper: 

Theorem 1.1. Assume that A^, n — >■ oo and that assumptions A\^and A\3[hold true. 
Assume moreover that there exists an integer p > 1 such that sup„E|Xj" |^P < oo and let 
{un) and (vn) be sequences of N x 1 deterministic vectors. Then, for every z e C — , 

E |< (Q„(z) - T„(z)) v^f' < ^$p(|z|)*, (^j-^ \\u,rP\\v4^^, (1.6) 

where (f>p and are nice polynomials depending on p but not on (u„) neither on (v„). 

Remark 1.1. Apart from providing the convergence speed 0{n~P), inequality (|1.6|) provides 
a fine control of the behavior of E\u* [Q — T)v\'^p when z is near the real axis. Such a control 
should be helpful for studying the behavior of the extreme eigenvalues of along the 

lines of [4] and [5]. 

Remark 1.2. Influence of the eigenvectors of AA* on the limiting behavior ofu*Qu. Consider 
a matrix S with no variance profile {D = In, D = In) and let T be given by (|1.4p . Matrix 
T writes in this case: 

Denote by V/S.V* the spectral decomposition of AA* , and by Ta : 

Ta.(-z(1 + .>+^) \ 

Obviously, T = VTaV* and by ThcorcmO u*Qu - u*VTaV*u . Clearly, the limiting 
behavior of u*Qu not only depends on the spectrum (matrix A) of AA* but also on its 
eigenvectors (matrix V). 
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Contents. In Scction[5J we describe two important motivations from electrical engineering. 
In Section [3J we set up the notations, state intermediate results among which Lemma 13.61 
which is the cornerstone of the paper. Loosely speaking, this lemma whose idea can be found 
in the work of Girko [11] states that quantities such as 

n 

u*Qiaia*QiU 

i=l 

are bounded. This control turns out to be central to take into account Assumption A-[21 An 
intermediate deterministic matrix i?„ is introduced and the proof of Theorem 1 1 . 1 1 is outlined. 
Basically, the quantity of interest u*{Q — T)v is split into three parts: 

u*{Q - T)v = u*{Q - EQ)u + w*(EQ - R)v + u*{R - T)v , 

each being studied separately. 

In SectionlU the proof of estimate of u*{Q—^Q)v is established, based on a decomposition 
of Q — EQ as a sum of martingale increments. Section [5] is devoted to the proof of estimate 
of ii*(E(5 — R)v\ and Section [6l to the proof of estimate of u*{R — T)v. 

Acknowledgment. This work was partially supported by the Agence Nationale de la 
Recherche (France), project SESAME n°ANR-07-MDCO-012-01. 



2. Two APPLICATIONS TO ELECTRICAL ENGINEERING 

Apart from the technical motivations already mentionned in the introduction. Theorem 
ll.ll has further applications in electrical engineering. In this section, we present an applica- 
tion to Multiple Input Multiple Output (MIMO) wireless communication systems, and an 
application to statistical signal processing. 

2.1. Optimal precoder in MIMO systems. A bi-correlated MIMO wireless Ricean chan- 
nel is a iV X 71 random matrix i7„ given by 

Hn = Bn + RI/'^—^^RI/'^ , 

where _B„ is a deterministic matrix, 1^ is a standard complex Gaussian matrix, and where 
Rn and i?„ represent deterministic positive N x N and n x n matrices. An important 
related question is the determination of a precoder maximizing the so-called capacity after 
mininum mean square error detection (for more details on the application context, see [I]). 
Mathematically, this problem is equivalent to the evaluation of a deterministic N x N matrix 
Kn maximizing the function T,nmse{Kn) defined on the set of all complex valued N x N 
matrices by 

N 

Imrase{Kn) = E^ [log (/ + Kr,HnW^Kl)]]\ (2.1) 
J = l 

under the constraint ■^Tr(iir„if *) < a (a > 0). This optimization problem has no closed 
form solution and one must rely on numerical computations. However, direct numerical 
attempts such as optimization by steepest descent algorithms or Monte-Carlo simulations 
to evaluate Xmmse{Kn) before optimization, or any combination of these techniques, face 
major difficulties, among which: Hardly tractable expressions for XmmseiKn), and for its 
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first and second derivatives, computationally intensive algorithms when relying on Monte- 
Carlo simulations. 

If N and n are large enough, an alternative approach consists in deriving a large system 
approximation Tmmse(-f^n) of Immse(^n) which, hopefully, is simpler to optimize w.r.t. Kn. 
This idea has been successfully developed in [T], in the case where i?„ = 0, and in in a 
slightly different context, where the functional under consideration is the Shannon capacity 
Is{Kn) = Elogdet (/ + KnH„H*K*). 

In the remainder of this section, we consider the case where i?„ 7^ and briefly indicate 
how Theorem 11.11 comes into play. First remark that for every deterministic matrix Kn , the 
random matrix K„Hn writes: 

where Wn is standard Gaussian random matrix (notice that (KnRn.K*)^^^'^ KnRl/^ is uni- 
tary). 

Using the eigenvalue/eigenvector decomposition of matrices KnRnK* and R„ , the unitary 
invariance of the canonical equations (|1.3p . and Theorem 11.11 one can easily check that 
the diagonal entries of (/ + A"„i/„i7,*A'*)~^ have the same asymptotic behaviour (when 
(ri, N) 00) as those of the deterministic matrix Tn{Kn) defined by: 

Tn{Kn) = [{I + ~S{Kn)KnRnK*) + + S{Kn)Rny' B*J<; 

where S{Kn) and 5{Kn) are the (imique) positive solutions of the system: 

5{Kn) = iTrX„i?„A-* [(/ + J(A-„)if„i?„A'*)+A'„B„(/ + 5(A-„)i?„)-iB*X;J 

(5(i^„) = iTri?„ [{I + 5{Kn)Rn) + BIK;-^{I +~5{K,-,)KnRnK:y^K,,B,^'^ 

(2.2) 

From this, it appears that T„-mse{Kn) can be approximated by Immse{Kn) given by: 

N 

TrnmseiKn) = ^ log [(/ + 6(Kn)K„R„K*) + + (5(A'„)7?„)"'b: A',*, 

Although the values taken by function Ar„ I„imse{Kn) are defined through the implicit 
equations (|2.2p . the first and second derivatives of Tmmse are easy to compute, and the 
minimization of Tmmse instead of T,nmse certainly leads to a computationally attractive 
algorithm. 

A number of important related questions remain to be addressed, e.g. the accuracy of the 
approximation Xmmse(^n), its impact on the error on the optimum solution, the derivation 
of a more accurate approximation as in [T], the development of an efficient algorithm to 
compute the optimal AT*, etc.; however this already underlines promising applications of 
Theorem 11.11 in the context of wireless communication. 

2.2. Statistical signal processing applications. There are many important applications 
such as source localization using antenna arrays, communication channel estimation, detec- 
tion of signals corrupted by additive noise, etc. where the observations are stacked into a 
matrix E„ given by p. II) in which An is a non observable deterministic matrix modelling the 
information to be retrieved and where Y„ is due to an additive noise. It is therefore often 
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relevant to estimate certain functionals of matrix A„ from I]„ . In this seetion, we show how 
Theorem 11.11 is valuable and relevant in the eontext of subspace estimators when N and n 
are of the same order of magnitude. 

Subspace estimation. Assume that ^ < 1, _D„ = In and -D„ = /„ (white noise); assume 
also that matrix Rank(^„) = r < N where r may scale or not with the dimensions n and 
N. Denote by n„ the orthogonal projection on the kernel of matrix An- The subspace 
estimation methods are based on the estimation of quadratic forms u*n„u„ where (u„)n6N 
represents a sequence of unit norm deterministic A'^-dimensional vectors. 

If TV if fixed while n — >■ +00, it is well known that ||I]„I]* — (A„yl* + /)|| 0. Hence, if 
tin represents the orthogonal projection matrix on the eigenspace associated to the N — r 
smallest eigenvalues of SnS* , then |1II„ — n„|| — > and thus 

u*n„u„ - u*n„u„ > . (2.3) 

n— >oo, N fixed 

In order to model situations in which n and N are large and of the same order of magnitude, 
it is relevant to look for estimators consistent in the regime given by (|1.2[) . Unfortunately, 
(|2.3[) is no longer valid in this context. 



An estimator for large N,n. The starting point of the estimator proposed in [25], inspired 
by [21], is based on the observation that n„ writes: 

n„ = / {AnA:^-xi)-' dX , 

where C~ is a clockwise oriented contour enclosing but not the non-zero eigenvalues of 
An An- In the white noise case, matrix Tn{z) writes: 

T„(z) = (1 + dniz)) {AnA*„ - Wn{z)iy^ , 

where w„(z) is the function defined by w„(z) = z(l + (5„(z))(l + (5„(z)). It is shown in 
[25j that (under additional assumptions) such a contour is the image under Wn of the 
boundary dTZy of the rectangle TZy = {z = x + iv, x € [x-,x+], \v\ < y} for well-chosen x- 
and x+ . A simple change of variable argument therefore yields the following formula for n„ : 

^n=^ I {AnAl ~ Wn{z)I)-' w\^{z)dz = / Tn{z) dz . 

Hence, u* n„u„ is given by: 

ullinUn = ^ f KTniz)Un ^'^'fl dz . (2.4) 

2*7r Jen- 1 + ''n(z) 

Eq. (|2.4|) is particularly interesting because all the terms in the integrand admit consistent 

estimators: Quantities (5„(z) and Sn{z) can be estimated by Sn{z) = ■^Ti'{Qn(z)) and Sn{z) — 

^Tr{Qn{z)), w^{z) can be estimated by the derivative of Wn{z) = z(l + 5n{z)){l + 5„(z)); 
finally. Theorem 11.11 implies that m* (5„(z)u„ — u*T„(z)u„ — >■ for A^, n — > oo. A reasonnable 
estimator for n„ should therefore be 

n„ = / Q„(z) /"f^\ dz (2.5) 
2«7r7a7j- l + Sn{z) 

and it should be expected that u*n„u„ — u*H„7i„ — > for N,n ^ oo. 
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Remaining mathematical issues. The full definition of n„ requires to prove that none of the 
poles of the integrand of the r.h.s. of (|2.5p can be equal to X- or x^. Otherwise, the mere 
definition of n„ does not make sense. This problem has been solved in the Gaussian case in 
|25j . In the non Gaussian case, partial results concerning "no eigenvalue separation for the 
signal plus noise model" [2] together with Theorem 11.11 tend to indicate that the estimator 
w*n„w„ is also consistent. 



3. Notations, preliminary results and sketch of proof 

3.1. Notations. The indicator function of the set A will be denoted by 1_a{x), its cardinality 
by #v4. Denote by aA6 = inf (a, b) and by aV6 = sup(a, b). As usual, = {x G M : x > 0} 

and C+ = {z e C : Im(z) > 0}; similarly C" = {z € C : Im(z) < 0}; i = \/^; if z e C, 

_ . . V 

then z stands for its complex conjugate. Denote by — > the convergence in probability of 

T> 

random variables and by — > the convergence in distribution of probability measures. Denote 
by diag(ai; 1 < i < k) the kx k diagonal matrix whose diagonal entries are the a^'s. Element 
of matrix M will be either denoted rriij or [M]ij depending on the notational context, 
if M is a n X n square matrix, diag(M) = diag(mii; 1 < i < n). Denote by M"^ the matrix 
transpose of M, by M* its Hermitian adjoint, by Tr (M) its trace and det(M) its determinant 
(if M is square). We shall use Landau's notation: By a„ = 0{bn), it is meant that there 
exists a nice constant K such that |a„| < iC|6n| as N,n — !• oo. Recall that when dealing 
with vectors, || • || will refer to the Euclidean norm; in the case of matrices, || • \\ will refer to 
the spectral norm. 

Due to condition (jl.2|) . we can assume (without loss of generality) that there exist < 
£~ < £+ < oo such that 

V7V,n e N*, £-< — <£+. 

n 

We may drop occasionally subscripts and superscripts n for readability. 

Denote by Y the N xn matrix n^^^^D^^^XD^^^; by (77^), (oj), (xj) and (yj) the columns 
of matrices S, A, X and Y. Denote by E^, Aj and Yj, the matrices S, A and Y where 
column j has been removed. The associated resolvent is Qj{z) = (SjS* — zIm)'"^ ■ Denote 
by Ej the conditional expectation with respect to the cr-field Fj generated by the vectors 
(y^i 1 < ^ £ i)- By convention. Eg = E. Denote by Ej,^ the conditional expectation with 
respect to the cr-field generated by the vectors (y^, £ ^ j). 



3.2. Classical and useful results. We remind here classical identities of constant use in 
the sequel. The first one expresses the diagonal elements of the co-resolvent; the other ones 
are based on low-rank perturbations of inverses (see for instance [TH Sec. 0.7.4]). 
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Diagonal elements of the co-resolvent; rank- one perturbation of the resolvent. 

Q(z)r]ir]*Q(z) 

QA^) = J'^ . (3-3) 



1 - 77;Q^: 



'J 



A useful consequence of 

Recall that ^2 = dist(2;, Considering the eigenvalues of Q{z) immediately yields 
I1Q(^)|| < • Taking into account the fact that 

and 



z(l + n-idjTr + a*Qjaj) ^(1 + 

are Stieltjes transforms of probability measures over R+, and based on standard properties 
of Stieltjes transforms (see for instance [15l Proposition 2.2]), we readily obtain the following 
estimates: 

^ < M and i r < M Vz e C - M+ . (3.6) 



l + ^Tri:iQj+a*Ojaj- 



The following lemma describes the behavior of quadratic forms based on random vectors 
(see for instance [H Lemma 2.7]). 

Lemma 3.1. Let x = (xi, • • • , a;„) 6e a n x 1 vector where the Xi 's are centered i.i.d. complex 
random variables with unit variance; consider p > 2 and assume that Ejxip^ < oo. Let 
M = (rjiij) be a n X n complex matrix independent of x. Then there exists a constant I'Cp 
such that 

E \x*Mx - Tr Mf < Kp (Tr MM*Y^^ . 
Let u e C" 6e deterministic, then E|a;*u|P = 0{\\u\\p). Moreover, ¥.\\x\\p = 0{nP/^) . 

Note by 7? = diag((ij ; 1 < i < N) and D = diag(c?i ; 1 < i < n). Gathering the previous 
estimates yields the following useful corollary: 

Corollary 3.2. Let z e C — IR+, and let p >2. Denote by Aj the quantity: 

A J = VjQjVj - ""TrDQj - a*Qjaj . 

Then 
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Theorem 3.3 (Burkholdcr inequality). Let (Xj.) be a complex martingale difference sequence 
with respect to the filtration {J-k)- For every p > 1, there exists Kp such that: 



E 



fc=i 



2p 



\k=l / fe=l 



A result on holomorphic functions: 

Lemma 3.4 (Part of Schwarz's lemma, Th.12.2 in }22j). Let f be an holomorphic function 
on the open unit disc U such that /(O) = and sup^gj/ \fiz)\ ^ 1- Then \f{z)\ < \z\ for 
every z €z U . 

Rules about nice polynomials and nice constants. Some very simple rules of calculus related 
to nice polynomials will be particularly helpful in the sequel: 

If ($fe,l < k < K) and (^'fe,l < k < K) are nice polynomials, then there exist nice 
polynomials <f> and 5* such that: 

K 

J2^k{^)'^k{y)<H^)'^iy) for x,y>0. (3.7) 
fc=i 

Take for instance $(a;) = J2k=i ^k{x) and ^{x) = J2k=i '^kix). 

If <i>i and ^'i are nice polynomials, then there exist nice polynomials $ and ^E* such that: 



V$i(a;)*i(j/) < $(x)*(y) for x,y > 0. (3.8) 
Take for instance $ = 2^^(1 + $1) and = (1 + v]/]^) and note that: 

1.. . ^ r ^,T, r ^^ / + 



V$l(x)*i(y) < -(1 + $i(x)^'i(y)) < ^ ^^^(1 + • 

The values of nice constants or nice polynomials may change from line to line within the 
proofs, the constant or the polynomial remaining nice. 



3.3. Important estimates. 

Lemma 3.5. Assume that the setting of Theorem ] 1.1\ holds true. Let u be a deterministic 
complex N X I vector. Then, for every z € C — R^, the following estimates hold true: 



p 



E|^E,_i(«*Qa,a*Q*«) I < i^p^l^ , (3.9) 



J2K,-i{u*Qv,rj*Q*u) \ < xj-^ip^ , (3.10) 



where Kp and Kp are nice constants depending on p but not on \\u\\ 



Proof of Lemma 15751 is postponed to Appendix El 



ON BILINEAR FORMS BASED ON THE RESOLVENT OF RANDOM MATRICES 



11 



Lemma 3.6. Assume that the setting of Theorem \l.l\ holds true. Let u be a deterministic 
complex N X I vector. Then, for every z G C — IR+. the following estimates hold true: 



J2K{u*Q,a,a*Q*uY < ^\z\)^S \\u\\\ (3.11) 

E (^E,_i(^.*Q,a,a*Q» < ^\z\)^ fj-\ \\u\\'p, (3.12) 



where (f>, ^P, $ and are nice polynomials not depending on \\u\\. 
Proof of Lemma is postponed to Appendix [K\ 

In order to proceed, it is convenient to introduce the following intermediate quantities 

(z e C-IR+): 

an{z) = -Ti- DnEQn{z), an{z) = -Tr DnEQn{z), (3.13) 

71 71 

Rniz) = (^-z(/Ar+a(z)I?„)/w + A„(/„ + a(z)A.) ' A*„^ , (3.14) 

Rn{z) = (-z(/„ + a(z)AO+A* (-^A^+5(z)A0~'^n) ' • (3.15) 
A slight modification of the proof of FT?} Proposition 5.1-(3)] yields the following estimates: 

\\Rniz)\\<^, \\Rniz)\\<^ forzeC-M+. 

Oz Oz 

The same estimates hold true for ||T„(z)|| and j|T!„(z)||. 

3.4. Main steps of the proof. In order to prove Theorem 11.11 we split the quantity of 
interest u*{Q — T)u into three parts: 

u*{Q - T)v = u*{Q - EQ)v + u*{EQ - R)v + u*{R - T)v , 

and handle each term separately in the following propositions: 

Proposition 3.7. Assume that the setting of Theorem ] 1.1\ holds true. Let (w„) and (w„) he 
sequences of N x 1 deterministic vectors. Then, for every z G C — M+, 

2P <r J_d>._nrlUT/._ ( —\ IU,_.I|2P|U, l|2p 



E |< (Q„(z) ~ EQ„(z)) < — $p(|z|)*p — ||7^ 



where and are nice polynomials depending on p but not on (u„) nor on (f„). 

Proposition 13 . 71 is proved in Section 21 
Proposition 3.8. Assume that the setting of Theorem \l.l\ holds true. 

(i) Let (un) and (u„) be sequences of N x 1 deterministic vectors. Then, for every 

|<(EQ„(z)-i?„(z))z;„| < -^Hlzl)^ (^) \K\\ \\vn\\, 
where $ and "i! are nice polynomials, not depending on (un) nor on (w„). 
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(ii) Let be a N x N deterministic matrix. Then, for every z £ C 



-Tril/„EQ„(z) - -TrM,,R,,{z) 
n n 



5_ 

where $ and "i! are nice polynomials, not depending on Af„ . 

Proposition 13.81 - (i) is proved in Section [5| proof of Proposition 13.81 - (ii) is very similar and 
thus omitted. 

Proposition 3.9. Assume that the setting of Theorem \1.1\ holds true. Let (u„) and (v„) be 
sequences of N x I deterministic vectors. 

Then, for every z € C — M.~^ , 

|< (i?„(z) - T„(z)) v,,\ < U{\z\)^ (^j-^ h,J |lt;„||, 
where $ and 4* are nice polynomials, not depending on (u„) nor on (w„). 
Proposition 13.91 is proved in Section [6] 



Theorem 1 1.1 1 is then easily proved using these three propositions together with inequality 
\x + y + z\^P < Kpi\x\^P + |ypP + |z|2p) and (|X7)) . 



4. Proof of Proposition 13.71 

Recall the decomposition: 

u*{Q - T)v = u*{Q - EQ)v + u*{EQ - R)v + u*{R - T)v . 

In this section, we establish the estimate: 
E|^*(Q(z)-Eg(z))«|^''<-L$p(|z|)*p(^l^||ri|l2f|K,fP, Vz e C - M+ . (4.1) 

4.1. Reduction to unit vectors and quadratic forms. Using a polarization identity, 
it is sufficient in order to establish estimate (|4.ip for the bilinear form u*{Q — E(5)u to 
establish the related estimate for the quadratic form u*{Q — ^Q)u and for unit vectors ||w|j 
(just consider w/||w|| if necessary): 

EK(g(z)-EQ(z))^|2^<-i<i>,(|z|)vI/p(^l) . (4.2) 

4.2. Martingale diflference sequence and Burkholder inequality. We first express the 
difference u*{Q — EQ)u as the sum of martingale difference sequences: 

n n 

u*{Q-EQ)u = ^(E, -E,_i)(w*gw) = ^(E, -E,_i)(u*(g-Q,)w) 
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One can easily cheek that {{Ej — Ej_i)rj) is the sum of a martingale difference sequence 
with respect to the filtration {J-j,j < n); hence Burkholder's inequality yields: 

2p 

^(E,-E,_i)r, 

3 = 1 



e(^5^E,_i|(E, -E,_i)r,f j +^^E|(E, -E,_i)r,f^ 



(4.3) 



Recall the definition of Aj — 'i]jQj'>]j — n ^djTr DQj — a*Qjaj. In order to control the 
right-hand side of Burkholder's inequality, we write Tj as: 

u*Qjr]*r]jQjU u*Qjr]*r]jQjU 1 + ^Tr DQj + ajQ^aj 



i- J — — — — — — X 



where 



1 + 7i*QjV3 1 + V^QjVj 1 + ^Tr DQ, + a*Qjaj 

1 + V*QjV3 1 + ^Tr DQ, + a*Q,a, 



^^^^^^^QfhVlQft r,, = ^ . (4.4) 

1 + i Tr Dg^ + a* Q.j f + ^ Tr DQ, + a* Qj 



In the following proposition, we establish relevant estimates. 

Proposition 4.1. Assume that the setting of Theorem \1.1\ holds true. There exist nice 
polynomials (<f>i, 1 < i < 4) and {'^i^ 1 < i < 4) such that the following estimates hold true: 

E('f]E,_i|(E,-E,_i)ri,|2)'' < l-<^^{\z\)^>J^\ , (4.5) 



f^E|(E,-E,_i)ri,f^ < 4'^2(|z|)^2f^) , (4.6) 



p 



E(5^E,_i|(E,-E,_i)r2,f j < -i-ci>3(|z|)vi/3 (^J- ) , (4.7) 
f^E|(E,-E,_i)r2,f^ = ^M\z\)^i(j-)- (4.8) 



It is now clear that the proof of Proposition 13.71 directly follows from Burkholder's in- 
equality together with the estimates of Proposition 14.11 The rest of the section is devoted 
to the proof of Proposition l4.1l 
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4.3. Proof of Proposition \4A\ Estimates (|i3|) and P^B)) . We split as = 
X2j + X3j, where: 



X2j 



u*Qjyjy*QjU 



l + ^TTDQ,+a*Q,a, 
y*QjUu*Qjaj 



a*QjUu*Qjyj 



1 + ^Tr DQj + a*Qjaj 1 + ^Tr ZJQ^- + a*Qjaj 



X3j 



l + ^TiDQj+a*Qjaj 



Notice that (Ej — Ej_i)(x3j) = 0, hence X3j will play no further role in the sequel. As Qj 
is independent from column yj, we have: 



1 + ^Tr DQj+a*Qjaj 



(4.9) 



and 



E,_i|(E,-E,_i)(xi,)l' < %^ 



X Ej_i 



c*i:'i/2Q^.yy*Q^.£)i/22;^. _ TrDQjUu*Qj 



l + ^Ti:DQ,+a*Q,a, 



r2 



max I I ^ 



< 



-2 



2- ^ ii^j-1 
2 



E, 



< i^^^ X E,_i fTri?l/2g ^^^*Q .^l/2^1/2g, ^*Q,^l/2 



O 



l2 <5! 



(4.10) 



where (a) follows from Jensen's inequality, (b) from estimate p. 61) . and (c) from Lemma |3. II 
Thus 

P / |^|2p 



E(^E,_i|(E,-E,_i)(Xi,)r 



(4.11) 



We now turn to the contribution of X2j- Arguments similar as previously yield: 



E,_i |(E, - E,_i)(x2,)r = E,-i |E,X2jr < ix2,r 



< -E,_i 
n 



2 \z\- 



]D^/'^Q.jUu*Qjaj 



1 + ^Tr i?Q,+a*Q,a, 



]QjUu*Q.jD^/'^Xj 



< -^-EJ_l (^Ej,^.(x*i:>i/2gjWw*Q*L»i/2a;j) x u*Qjaja*Q*u 



+E.y^.(a-*i:'i/2Q*^^*Q^.^i/2^.^.) ^ u*Q*aja*Qju) , 

(4.12) 



^' |2|2 

< {u*Q*aja*Qju) +Ej_i {u*Qjaja*Q*u)) 
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Now, using Eq. p.l2p in Lemma [3.61 yields : 



(4.13) 



We now establish estimate (|4.6p . As previously, consider identity (|4.9p : take it this time 
to the power p. Using the same arguments as for (|4.10p . we obtain: 



Hence, gathering (|4.1ip and (|4.13p yields estimate (|4.5p . 



hence: 

EEl(E.-IE.-i)(Xi.)r^ = o{-J^) ■ (4-14) 

Similarly, using the same arguments as in (|4.12p . together with elementary manipulations, 
we obtain: 

E,_i |(E, - E,-i)(x2,)l'' < §^ {u*Q*a,a*Q,uY + {u*Q,a,a*Q*uY) . 

Due to the rough estimate (jA.ip . we obtain 

E |(E, - E,_i)(x2j)l''' < 4 j£ (e {u*Q*a,a*Q,uy + E {u*Qja,a*Q*uY) , 
which after summation, and the estimate obtained in Lemma |3.6[ yields: 



Ef]|(E,--E,._i)(X2,)|'^ < ' (4-15) 

where <!>' and are nice polynomials. Gathering (|4.14p and (|4.15p yields estimate 



4.4. Proof of Proposition llTTt Estimates and We split T2j as T2j = Xij 

X2j + X3j, where: 



Xij = X 
X2j - Aj- X 
X3j = Aj- X 



u*Qjaja*QjU 



(1 + 77*Q,?/j)(l + ^Tr i?Q, + a*Q,a,) 



(l + 77*Q,?/,)(l + ^Tri^Q,H 








(1 + ^*Q.?7,)(l + ^Tr Z?Q,H 




u*Qjyja*QjU + u*Qjajy 


jQjU 



16 



HACHEM ET AL. 



Consider first: 



E,_i|(E, -E,_i)(xi,)r < 2E,-i|xi,f 



< 4 Ej„i u*Qjaja*QjU lyjQjUj ~ n djTrDQj 



-Ej_i \u*Qjaja*QjU (y*Qjaj + a*Qjyj)\ 



(b) K\z\' 



< 



< 



i2 6l 



K\;. 



n S 



K\zl 



Ej_i u*Qjaja*Q*u Ey.{x*D^^^Q*aja*QjD^^^Xj) 



W K\z\^ 



Ej-i{u*Qjaja*Q*u) , 



where (a) follows from p.6p . (6) from the fact that \u*Qjaja*Qju\ < KS^ ^ and \u*Qjaja*Q*u\ 
KS^^^ and (c) from Lemma |3. II From this and Lemma |3.6[ we deduce that: 



< 



E |^E,_i|(E,-E,_i)(Xy)r 



(4.16) 



Consider now: 



(a) 



E,_i|(E,-E,_i)(x2,)r < 2E,_i|x2jf < ^E,_i |A,r < 



^3 ' 



where (a) follows from the triangle and Jensen's inequality, (b) from p.6p and (c) from 

Cauchy-Schwarz inequality, Lemma |3. II and Corollary 13.21 

Hence, 



e|^E,_i|(E,-E,_i)(x2,)|' 



Similarly, one can prove that: 



E |^E,_i|(E,-E,_i)(x3j)r 



O 



Gathering the previous results yields the bound: 



E(^E,_i|(E,-E,_i)(r2,)|'| < 
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Wc now evaluate the second part of Burkholdcr's inequality (and may re-use notations $ 
and 4* for different polynomials). 

^E|(E,.-E,_i)(Xi,)l'" < -^Y.^{u*Q,a,a*Q*uf''Ey^\^,\^^ 
■j=i ^ i=i 



< 



^ ^ E {u*Qjaja*Q*u) {u*Q.jaja*Q*uy 



z j = l 



2 J = l 

< J_$(U|)^fJ_ 



where (a) follows from Corollary 13.21 and the last estimate, from Lemma 13.61 Similar com- 
putations yield: 

^E|(E,-E,_i)(X2,)|'^ < ^,^\\z\w(^\ , 

" 1/1 

^E|(E,-E,_i)(X3,)l'" < 

j=i ^ ^ 

the first of these inequalities requiring the assumption sup„ EjJ'Tj" < cx) in the statement 
of Theorem ll.il Gathering these three results yields: 

X:E|(E,-E,_0(r2,)l^^ < ^nA)^[^ , 

j=l \ 

and Proposition 14. II is proved. 



5. Proof of Proposition 13.81 

Recall the decomposition: 

u*{Q - T)v ^u*{Q~ EQ)v + u*{EQ - R)v + u*{R~ T)v . 
In this section, we establish the estimate: 

1^.* (EQ(z) - R{z))v\ < ^$(kl)* (^) \\u\\ \\v\\ , 

The argument referred to in Section (|4.ip still holds true here; therefore it is sufficient to 
establish, for z G C — R+ and for a unit vector u: 

\u* (EQ(z) - R{z)) u\ < -^^(kl)* (jA , (5.1) 
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r ~ 1 

Recalling that R = —z{I + aD) + A{I + aD)^^ A* , the resolvent identity yields: 



u*{R - Q)u = u*R{Q-^ - R-^)Qu , 

= u*R (SS* - A{I + aOy^A* ) Qu + zauRDQu 



= u*R \ V ?7,?7* - V Qu + zau*RDQu 

(a) ■J^ u*Rrijri*QjU u*Raja*QjU 



j=i \- ' 'ij^riuy- ' — j=i 



A 



where (a) follows from (|3.2p and p.Sp . together with the mere definition of a. 

As usual, we now write i]j = + aj , group the terms that compensate one another and 
split Zj accordingly: 

where 

y*QjUu* Ri/i dn I 1 \ 



^2, = 



^3 J 



{adj - y*Qjyj)u*Raja*QjU 

y*QjUa*Qjyj X u*Raj 

u* Ryja*QjU + u* Rajy*QjU 

y* Qj ajU* Raj a* Qj u + a* Qj yjU*Raj a* Qj u 

u*Raja*Qjajy*QjU + u*Raja*Qjyja*Qj 
{l+V*QjVj){^ + adj) 



Now, the estimate (|5.ip immediately follows from similar estimates for the terms E 
1 < £ < A. The rest of the section is devoted to establish such estimates. 
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5.1. Convergence to zero of J^j^^ij. We have 



KZij ^ E 



E 



y*QjUu*Ryj \ _ dj_^ ( 1 
^ 1 + VjQ]Vj J n \l + ri*QjT]j 

/y*Q^uu*Ryj\ f u*RDQ,jU^ 



E{u*RDQu) 



n 



l{u*RDQju) 



A 



Xlj + X2j + X3i ■ 



¥.{uRD{Qj - Q)u) 



We first handle . RecahthatAj = rjjQjr/j—n ^djTr DQj~a*Qjaj. Since E,y .{y*QjUU* Ryj) 
djU^^u* RDQjU, we get: 



Xlj 



E 



E 



E 



y*QjUu*Ryj\ d^ f u*RDQjU 
1 1 



A, 



1 + Tl*,QjTlj 1 + iTr DQj + a*Q,aj 

y*QjUu*Ryj - ^{u*RDQju) 
(1 + V*Qjm){^ + ^TrI?Q, + a*Q,a,) 



y*QjUu*Ryj '-{u*RDQju) 



Hence, 



1X1,1 < 



E 



< 



y*QjUu*Ryj - -^{u*RDQ^u) 



= O 



1/2 



8l V^S, nS'l ^ KnV^S 

Summing over j yields the estimate Ixijl = 



-5^ 



We now handle X2j- Using the inequality cov(Xy) < ■\/var(X)var(y), we get: 



, , K\z\ 



\u*RD{Qj - EQj)u\ 



Hence, applying Proposition 13.71 to \u* RD{Qj — MQj)u\'^ and summing over j yields the 
estimate J2j 1x2,1 = ?i-i/2$(|z|)*(^7'). 
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Let US now handle the term xsj- Using the decomposition of Qj — Q, Schwarz incquahty 
and the fact that \/ab < 2^^ {a + b) yields 



1 



K\z? 



E{u*RD{Qj - Q)u) 



< {E\u*RDQ,,y\'+E\r^*Q,u\^) 



(5.2) 



Now, as: 



E\u*RDQjr]j\^ = Eu* RDQjyjy*Q*DR*u + Eu* RDQjaja*Q*DR*u , 
E\ri*Qju\'^ = Eu*Q*yjy*QjU + Eu*Q*aja*QjU , 



it remains to sum over j and to apply Lemma [3.6l to get the estimate | xsj \ — n ^^(\z\)'ii(S^^). 
Gathering the partial estimates yields: 



eY^z 



< 



H\z\)^is-') 



(5.3) 



5.2. Convergence to zero of J^j^^^j- Recall that 



Z2j - 



[adj - y*Qjyj)u*Ra.ja*QjU 



We have: 



\EZ2j\ < 



^ ^\u*Raj\E {adj - y*Qjyj)a*QjU 



< ^\u*Raj\^E\a*Qju\^\lE adj - y*Q,yj 



< 



/ u*Raja*Ru + Eu*Qjaja*QjU 



E 



Ydj - y*Qjyj , (5.4) 



where (a) follows from p.6p . In order to estimate the remaining square root, we decompose 
the difference as: 

adj - y*Qjyj = ^Tri?(EQ - Q) + DiQ - Qj) + ^Ti DQ^ - y*Q,yj . 

Hence, 

E\adj - y*Qjyjf 



<K I -^E|TrD(EQ-Q)|^ + ;;^E|TrD(Q-0j-)|2+E 



'k^,DQ,~y*Q,y, 



Writing E|n-iTr£i(Q-EQ)|2 < 1+ supj<„ E|e*D(g-EQ)ejf where Cj represents canon- 
ical vector number j and using the result of Section HI the first term of the right hand side is 
of order n~^(f>(|z|)5'(^J^). The second term is of order {nSz)"'^ (minor modification of [2H 
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Lemma 2.6] to encompass the case Re(z) < 0). Finally, the third term is of order n ^5^ ^ 
by Lemma [3TTJ Collecting these results, we obtain: 

1/2 



K 



$2*2 



1/2 



(a) K (b) K 

where the $'s are nice polynomials with argument \z\ and the ^''s are nice polynomials with 
argument \S~ |, and where (a) follows from p.7p and (b) from (|3.8p . It remains to plug this 
estimate into (|5.4p . to sum over j and to use Assumption [5] together with Lemma 13.61 to 
obtain: 



< 



K\zl 



< _$'(|„.|)vl/'(5;i) 



(5.5) 



5.3. Convergence to zero of J^j^^sj- Recall that 



3j 



y*Qjua*Qjyj x u*Raj 
{l + j^*Q,il,){l + adj) 



We have: 

E|Z3,| < ^-^\u*Ra,\xny*Q,ua*Q,y,\ < \u*Ra, \ ^E\y*Q,u\m\a*Q,y,\^ 

^ Z ^ Z 

< -^\u Raj\ , 

n 

where (a) follows from (|3.6I) . and (b) from Lemma |3. II Hence, 



< —^y"\u*Ra,\ 



n S 



\ 



'^u*Raja*R* 



= O 



3 = 1 



IS' 



(5.6) 



5.4. Convergence to zero of J^j^^-ij- Write Z4j 



{l + 'n*Q,r,j){\ + adj) 



with 



Wj^j = (1 + adj){u* Ryja*QjU + u* Rajy*Qju) 



y*QjajU* Raja*QjU — a*QjyjU*Raja*Qj 



+ u* Raja*Qjajy*QjU + u* Raja*Qjyja*QjU 
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Write 



A, 



1 + V*Q3V3 1 + ^TrZ?Q, + a*Q,a, (1 + v*QjVj)i^ + ^Tri^Qj + a*0,a,) ' 
Plugging this identity into and taking into account the fact that Ky. W^j = 0, we obtain: 



E 



(1 + adj)il + V^QjVj){^ + TfTr-DQj + ^QjOj 



K \z\'' 



< ^^E|A,iYE|i^4,p < -^^^/m^,\ 



Hence, 



< 



We therefore estimate J^j I^|W^4jP- First, write: 



E|W4jf < -^U + y] ma*Qju\^u*RDR*u+\u*Ruj\^E{u*Q*DQju)) 



K 



Now, summing over j yields: 



$^£1^-4,12 < ^ ( ^E(w*g*a,a;Q,«) | ("l + i-) 1 



< i4.(| = |)(I>(5ji) . 



Plugging this into (|5.7p yields the estimate 



Jn 



(5.7) 



(5.8) 



5.5. End of proof. Recall that: 
\u*{R-EQ)u\ < 



n 

E^Zi, 


+ 


n 


+ 


71 


+ 


n 

E^Z4, 



It remains to gather estimates (|5.3p . (|5.5p . (|5.6p and (|5.8p to get the desired estimate 

\u*iR^EQ)u\ < -^mz\)^{S;') . 
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6. Proof of Proposition 13.91 

Recall the deeomposition: 

u*{Q - T)v = u*{Q - EQ)u + u*(EQ - R)v + u*{R - T)v . 
As mentioned in Section HTTl it is sufBcient to establish the estimate: 

\u*{R{z)~T{z))u\<^nz\)^l^£^ , (6.1) 

for z G C — in the case where u has norm one. 

6.1. The estimate for u*{R-T)u. Recall the definitions of 5, ~5 a, a (|3A3l) and R, R 
(|3. 14113. 15| . Using twice the resolvent identity yields: 

u* {R -T)u^{a- d)Ki + {a-d)K2 , (6.2) 

where 

Ki = zu* RDTu 

K2^u*RA{I + aD)-^D{I + SD)-^A*Tu ' 
The following bounds are straightforward: 

|ki|<^%^ and |k2| < M^J^^ X ||(/ + a^)-i|| X ||(/ + 5^)-i|| . 

It remains to control the spectral norms of (/+ al))"^ and (/ + SD)~^. Recall that a is the 
Stieltjes transform of a positive measure with support included in R"*". This in particular 
implies that Im(zQ;) > for z e C+. One can check that 

Tj(z) = i ^ 

— z(l + ctdj) 

is analytic and satisfies Im(Tj) > and Im(zTj) > on C''^ and that limy^oo{—iy'^j{iy)) = 
1. As a consequence, Tj is the Stieltjes transform of a probability measure with support 
included in (see e.g. [151 Prop. 2.2(2)]). In particular, 

|T,(z)|<^ for l<j<n, 

which readily implies that ||(/ + < The same argument applies for ||(/ + 

SD)~^\\. Finally, 

|K2| < 74 ■ 

In view of the estimates obtained for ki and K2, it is sufficient, in order to establish (|6.ip . 
to obtain estimates for a — S and a — S. Assume that the following estimate holds true: 

VzeC-M+, max(|a-(5|,|a-5|) < i$(|z|)* , (6.3) 



where $ and are nice polynomials. Then, plugging (|6.3p into (|6.2p immediately yields the 
desired result (16.11). 



The rest of the section is devoted to establish 
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6.2. Auxiliary estimates over (a— (5) and ((5— (5). Writing a = n ^TrDR+n ^TiD{EQ— 
R) and(S = n'^TrDT, the difference a - (5 expresses as n-^Tr D{R~T) + n-^Ti- D{EQ - R). 
Now using the resolvent identity R — T — ~R(R~^ — T~^)T and performing the same com- 
putation for the tilded quantities yields the following system of equations: 



a — 5 \ ^ 1 ct — 5 \ ( e 



. , , = Co ( : ?) + (-) where Co = ( I"" ) , (6.4) 

a- 5 J \ a - 6 J \ e J \ zvq uq J 

the coefficients being defined as: 

uo = ^TtD^/^RA{I + al)y^b{l + 5b)~^A*TD^/'^ 
uo = lTrD^^^RA*{I + aD)-^D{I + 6D)-^AfD^/^ , . 

vo = '^TrDRDT ' ^ ' 

^ n 

= ^TyDRDT 

^ n 

and the quantities e and i being given by: 

e = -TyD{EQ-R) and e = -T^x b{¥.Q - R) . (6.6) 
n n 

The general idea, in order to transfer the estimates over e and i (as provided in Proposition 
13.81 -fii) ), to a — (5 and a — (5, is to obtain an estimate over 1/ det(/ — Co), and then to solve 
the system (|6.4I) . 

Lower bound for det(/ — Cq). The mere definition of / — Co yields 

|det(/-Co)| = \{l-uo){l-uo)- z^vqVo] 

> (l-|uo|)x(l-|uo|)-kn«o|x|f}o| 

In order to control the quantities uo,uq,vq and vq, we shall use the following inequality: 

\Ti AB*\ < {Tr AA*)^^^ X {TtBB*)^^^ , (6.7) 

together with the following quantities: 

^Tr DTA{I + 5*b)-^b{I + 5b)-^A*T* 
lTrbfA*{I + SD)-^D{I + ~5*D)~^Af* 
'^TrDTDT* 

n 

^TyDTDT* 



Ul 
Ul 

Vi 



and 



U2 = ^TrDRA{I + a*by^b{I + ab)-^A*R* 
U2 = lTrbRA*{I + aD)-^D{I + &*D)-^AR* 
\tvDRDR* 

"^TibkbR* 



(6.8) 



V2 

Using (HJI) together with identity (/ + 6b)-^A*T = fA*{I + SD)''^ (and similar ones for 
related quantities), we obtain: 

l^^ol < iuiu^y^' , l^iol < iuiU2)'/^ , li^ol < iviV2)'/^ , \vo\ < ivii2)'/^ , 

hence the lower bound: 

I det(/ - Co)| > (1 - {uiU2y^'){l - {u,U2f'^) - \zWv,V2iii2f/^ . (6.9) 

Notice that it is not proved yet that the right hand side of the previous inequality is non- 
negative. 

In order to handle estimate (16. 9p . we shall rely on the following proposition. 
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Proposition 6.1. Consider the nonnegative real numbers Xi,yi, Si,ti (i — 1,2). Assume 
that: 

Xi<l, yi<l and (1 - x^){l - yi) - Siti > for i = 1, 2. 

Then: 



(1 - ^XiX2) (1 - y/mh) - y/siS2tit2 



> v/(l - - yi) - Siti^{l - X2){1 - 2/2) - S2t2 . 

Proof. If a > c ( > 0) and 6 > d ( > 0), then: 

\/ab — \fcd > \/a — c\Jh — d . 

To prove this, simply take the difference of the squares. Applying once this inequality yields 
1 - y/xTx^ > ^/(T^^xi)(l^X2) , hence: 



(1 - y^XiX2) (1 - y/mm) - VsiS2hh > V(l ^ 2;i)(l - X2)il - yi)(l - 1/2) - VsiS2tit2 

Applying again the first inequality yields then the desired result. □ 

Our goal is to apply Proposition l6.1l to (j6.9|) . The main idea, in order to fulfill assumptions 
of Proposition 16.11 (at least on some portions of C — K'*'), is to consider the quantities of 
interest, i.e. Ui,Ui,Vi,Vi (i — 1,2) as coefficients of linear systems whose determinants are 
the desired quantities (1 — Ui)(l — Hi) — \z\'^ViVi. 

Consider the following matrices: 

C,;(z)= ( , V 1 = 1,2. 

The following proposition holds true: 
Proposition 6.2. Assume that 2: e C — M+. Then: 

(i) The following holds true: 1 — ui{z) > and 1 — ui{z) > 0. Moreover, there exists 
positive constants K, 77 such that: 

det(/-Ci(z)) > Kj-^^^^. 

(ii) There exist nice polynomials $ and and a set 

f« = |^eC+, l$(|z|)vi/ (^J-^ < 1/2 

such that for every z £ £n, 1 ^ ^2(2) > 0, 1 — U2{z) > 0, and 

5« 



det(/- C2) > K- 



(ry2 + |z|2)4 
where K, rj are positive constants. 



Proof of Proposition 16 . 21 is postponed to Appendix iBl 

We are now in position to establish the following estimate: 

yzeSn, max(|a-<5|,|5-^|) < ^$(|z|)* (^j . (6.10) 
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Assume z €z £n- Thanks to Proposition 16.21 assumptions of Proposition 16.11 arc fulfilled by 
Ui,Ui,Vi and Vi, and (|6.9p yields: 



dct(J-Co) > v/dct(/-Ci)Vdet(/-C2) > K / , (6.11) 

where K, rj are nice constants. 

Solving now the system ()6.4p . we obtain: 

a— 5 = {dct{I — Cq))^^ {{1 — ua)e + zvoi) 



a — 5 = (det(/ — Co)) ^ ((1 — uo)e + zvoe) 

It remains to use ()6.1ip . Proposition 13.81 - (ii). and obvious bounds over uo,uo,vo and vo to 
conclude and obtain (|6.10p . 

We turn out to the case where z G C — — £„, and rely on the same argument as in 
Haagerup and Thorbjornsen [T2] (see also [7|). In this case, 

U{\z\)^is:')>^. 

As \a-S\ = \n-^TrD{EQ - T)\ < 2£+d^^^S~\ we obtain: 

r,f+. 2$(|z|)^'fj- 



a similar estimate holds for a — 5 for z ^ En- Gathering the cases where z G £n and z ^ En 
yields 



Appendix A. Remaining proofs for Section O 

Proof of Lemma 13.51 Note that it is sufficient to establish the result for a vector u with 
norm one (which is assumed in the sequel). The general result follows by considering 

We proceed by induction over -p. Let p = I and consider: 

n 

< E^Ej.iu*Qaja*Q*u = Eu*QAA*Q*u < a^^^^E\\Q\\^ . 



As HQ II < ^, we obtain the desired bound. 
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E 



Now, write 

n 

'^Ej^i{u*Qaja*Q*i 

3 = i 



= E^E^^-iiu*Qaj^a*^Q*u) ■ ■ ■Ej^^i{u*Qaj^a*^Q*u) 

ji,--- dp 

< p! ^ E E.i^^i{u*Qaj^a*^Q*u) ■ ■ ■Ej^^i{u*Qaj^a*^Q*u) 
ji<---<ip 

p-i 

J-'jp - 1 measurable 
n p—l 



h<---<ip-i 3p=jp~i 

2 



J t '-''max 



E 



-1 



where (a) follows from the fact that 



0p=3p-l ip = ^ ^ 

It remains to plug the induction assumption to conclude. Hence p.9p is established. 

In order to establish p.lOp , one may use the same arguments as previously together with 
the identity QSS* = / + zQ, which yields the factor |z|^ in estimate (|3.10p . 

Proof of Lemma 13.61 We prove the lemma in the case where j|u|| = 1, the general result 
readily follows by considering 

Write u*Qjaja*Q*u = Xij + X2j + Xsj + Xij with: 

Xij = u*{Qj - Q)aja*{Qj - Q)*u 

X2j = u*Qaja*Q*u 

X3j = u*{Qj - Q)aja*Q*u 

Xij = u*Qaja*{Qj - Q)*u 

Hence, 



Notice that: 

E|x3,f < I (Ex?, + EX2j) and E|x4,f < ^ (Ex?, + Ex^,) • 
Note that using the facts that aja* < AA* and 77,77* < SS* together with the identity 



QSS* = / + zQ yield the rough but useful estimates: 

u*Qaja*Q*u = O (<5;^) and u*Q'qj'q*Q*u = O ( ^ 



(A.l) 
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Wc first begin by the contribution of ^ Exii 



^X2j = '^^u*Qaja*Q*u X u*Qaja*Q*u , 

n 

< '^u*Qaja*Q*u x u*QAA*Q*u , 

< {u*QAA*Q*uf = 0{S-'^) 

< $2(kl)*2 ^ 



Similarly, 



We now turn to the contribution of ^Xij ■ Using the decompositions 
p.4p . xii writes: 



1 + ?7*Qj77j 



1 - VjQV] 



X \u*Qrijri*Qaja*Q*r]jr]*Q*u\ 



1 - r]*Qr]j 



We first prove that 



a*Q*7jj7j*Qai 



1 - r]*Qr]j 



O 



In fact: 



1 - ?7*Q??j 



< 



(a) 



ajQ*VjVjQ*aj 



1 - 



< \a*{Qj - Q)*aj\ + 2|Im(z)||a*(Q, - Q)Qa 



where we use the fact that Q — Q* ~ 2ilm{z)Q*Q to obtain (a). Now, 



^TtDQj +a*Qja, 
n J J J 
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Since \n~^djTT DQj +a*Qjaj\ = 0(Sl^^), wc obtain: 

'"''^ ■-o(\^))xJ2^{u*Qr^,rf^Q*uY 



o 



(a) 



|z|4 



< $l(l2l)*l 



^1° 



where (a) follows from (|A.3[) and (lA.ll) and (6), from Corollary 13.21 
It remains to gather the contributions of XijTX2j,X3j ^^nd Xij to get: 

J2K{u*Q,a,a*QjuY < 2<i>,{\z\)^, ( j-) +2<!>2{\z\m ( j-) < $(kl) * f 



3=1 ^ 

where (a) follows from p.7p . Eq. p.lip is proved. 
In order to prove (|3.12p . first note that: 

< I E 



E 



p 


n 


p 




n 




+ E 






KE 







Hence, it remains to evaluate the contributions of each term. Using decomposition (|A.4I) 
together with the estimate (jA.Sp . we obtain: 



E 



- O ( ^ j X e(^ ^E,_i|l + ^]Q,^,\ X u*Q^,^*Q*u 



Using (|A.6|) together with (|3TT0)) yields 



E 



|2p 



o 



X E 



^E,_i(|A,| xw*g7;,77*g* 



Combining standard inequalities (Cauchy-Schwarz, \Yl,^a^hj\ < ^f)^^^' ^'^'^ 

Cauchy-Schwarz again), we obtain: 

n 
3 = 1 



< 



e( J2^3~ii^*Qv3VjQ 

3 = 1 



*U? I X 



e( ^E,_i|A, 

3 = 1 



1/2 



(a) 
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where (a) follows from (|A.1|) . Corollary and ([5TTU|) . Finally, 

P /l,|2p\ /U|2p\ /U|2p 



E 



Eq. f|3.9p directly yields the estimate: 



O ( 4^ ] +0(^4^] +0 ( ^ ) < $i(|z|)*i(^;^) . (A.7) 



o 



Finally, 



1/2 



E 



< E 



i=i j=i 



< $3(|^|)*3 



A corresponding inequality exists for E| '^^j-iXijl^'- obtain: 



E 



^Ej„iX4j 



< $4(k|)*4 



(A.8) 



(A.9) 



(A.IO) 



Gathering (|Aj|, (KM, (SM and ([AlOl) . we end up with p.l2p . and Lemma [3^ is proved. 



Appendix B. Remaining proofs for Section [6] 

Proof of Proposition [HI]- (i). RecaU that 5 = j^Tr DT and 6 = ^TrDT. We consider 
first the case where z £ C+ U C^. We have 

Im((5) = —Tt DT(T-* -T-^)T* and Im(zJ) = — Tri)(zf) \(zf)~* - (zf)"!] (zf)* . 
2in 2in L J 

Developing the previous identities, we end up with the system: 
where 



(B.l) 



wiiz) = ^TrDTT* (> 0) 

= ^TrDfA*{I + 6D)-^{I + S*D)-^Af* (> 0) ' 

By developing the first equation of this system, and by recalling that S{z) is the Stieltjes 
transform of a positive measure fin with support included in R+ , we obtain 

lm(z) lm(z5) lm(z) 

Replacing {lm{5),lm{zS)) with (Im((5), Im(z(5)) and repeating the same argument, we obtain 

1 . - Ini(z) , Im{zS) ^ Im(z) 

lm{S) lm{S) ~ lm{S) ~ 

By continuity of Ml (z) and{2i(2;) at any point of the open real negative axis, we have > 

and 1 — Ml > for any z G C — M+. The first two inequalities in the statement of Proposition 
(i) are proven. 
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By applying Cramer's rule ([ini Sec. 0.8.3]) where the first column of / — Ci is replaced with 
the right hand member of (jB.ip . we obtain 

„ s , N Imfz) Im(z) , , Im(z) Im(z) Im(z) ,„ ^, 

Using the fact that the positive measure is supported by IR+ and has a total mass ??^^Tr Z?, 
we have 



Im(,5) f 1 1 1 T n <- ^^'^'"ax , „ ^ Im(^) d, 



max 



Im(z) J |t-z|2^"^ Sin - Si ' " Im(z) " St 

' (B.3) 

In order to find a lower bound on wi and we begin by finding a lower bound on \6\. 
A computation similar to [151 Lemma C.l] shows that the sequence of measures {fin) is tight. 
Hence there exists r/ > such that: 

Mn[0,.y]>i-Tri?>^l^. 
2 n Z 

We have 

" '4(?/2 + |2|2) 

(B.4) 

Furthermore, when Re(z) < 0, we have 



|<S| > Re(<S) = J L-^^,^(dt) > -Re(z) j > ™ Re(z) 

|<5|>^ 



4(r;2 + |z|2) 



which results in 

^ "min 

^4(772 + |z|2) ■ 

We can now find a lower bound to wi: 

N N N N 

wi = -Ti-DTT* = -^d,"^\T,jf = -Tri?^K,^|ryf with k 



n n ^-^ ^-^ n ^-^ ^-^ ir 

i=i j=i 1=1 j=i 



> ^Tri?|^«Jf:|T,|A j > iTri?(^^.,|T,|^ >iTri? 



TV 



iTri^ - 16 £+d„ax(r?2 + |z| 



2^2 



where (a) follows by convexity. A similar computation yields wi > ((52dmin)^/(16 dmax('7^ + 
|z|2)2) where 77 is a positive constant. Grouping these estimates with those in (|B.3p and 



plugging them into (|B.2p . we obtain 



det(/- Ci) > 



*min / 



256 (£+d^axd„,ax)2(772 + \z\^)Hlf + \z\^)^ 

> K- 



8l 



(max(7/, 77)2 4- |z|2)4 

where if is a nice constant. The same bound holds for z £ (—00, 0) by continuity of det(/ 
Ci(z)) at any point of the open real negative axis. 
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Proof of Proposition [6T21 - (ii) . Recall that 

1 



e„ = -TrD(EQ - R) 
n 



We first establish useful estimates. 

Lemma B.l. There exists nice polynomials $ and such that: 



Im(£„(z)) 



Im(z) 



f 



< -^(Izl)* 



f 



and 



Im(ze„(z)) 



Im(z) 



n 



for z G C - 



(u), the sequence of functions (e„) 



Proof. We prove the first inequality. By Proposition 
satisfies over C — M+ 

\en{z)\ < -$(|z|)f 
n 

where $ and 5* are nice polynomials. Let TZ be the region of the complex plane defined as 
11 = {z : Re(z) < 0, | Im(z)| < - Re(z)/2}. U z & C-R+ -TZ, then | Im(z)| > <5Jv^, 
therefore |Ime(z)/lmz| < n~^\/5S~^^{\z\)'${S~^) and the result is proven. Assume now 
that z € TZ. Li this case, z belongs to the open disc centered at Re(z) with radius 
-Re(z)/2. For any u € V^, we have |e(u)| < n~^<I>(|u|)*(|u|"i). Moreover, 



„ 6z Re(z) , , 
Vug 2?,, -^< ^<l"l< 



3Re(z) ^ 3|z| 
2 - ~ 



As <J>(a;) is increasing and 5'(1/.t) is decreasing in a: > 0, we obtain: 



1 



3|z| 



\eniu)\ < -$ l^-^j * for ueV, . 

The function e is holomorphic on . Consider the function: Applying Lemma 

^ £(|Re(z)/2|C + Re(z))-£(Rc(z)) 
sup„gpJeH-£(Re(z))| 

Let C ~ 12 Im(z)/ Rc(z), apply Lemma [3.41 and use (jB.5[) . This yields: 

2|Im(z)| \/5|Im(z)| 
|Re(z)| n \0z J 8z n 



(B.5) 
with 



|£(z)-£(Rc(z))|< 



where $ and 4* are nice polynomials. As Im(£(Re(z))) = 0, we obtain 



Im(£„(z)) 



Im(z) 



< 



£(z)-£(Re(z)) 



Im(z) 



<^$(|z|)vl/ 



This proves the first inequality. The second one can be proved similarly. 



□ 



We now tackle the proof of Proposition 16 . 21 - (ii) . following closely the line of the proof of 
Proposition ElKi). Recah that a = ^TrDEg, 5 = iTrI)EQ, e = iTrD(Eg - R), and 
£ = iTrI)(E(5 — R). We begin by establishing the lower bound on det(/ — C2). Assume 
that z G C+ U C~. Writing a = ^Ti DR + e and a = ^Ti DR + e and developing Im(Q;) 
and Im(za) with the help of the resolvent identity, we get the following system: 



lm(Q;) 
Im(za) 



Im(z) 



W2{z) 

X2{z) 



Lir(£) 
Lii(z£) 
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where W2{z) — ^TrDRR* and X2{z) > 0. Let W2 ~ n ^Ti DRR* . Using the same argu- 
ments as in the proof of Proposition 16. 2l -(i). we obtain 

^ ^ Im(z) ^ Im(za) ^ Im(e) ^ Im(z) ^ Im(£) 

lm(a) Im(a) Im(a) ~ Im(a) Im^a) ' 

2 - ~ Im(z) ^ _ Im(za) ^ Im(£) ^ _ Im(z) ^ Im(£) 

lm{a) Im(a) Im(Q!) ~ Im(Q!) lm{a) ' 

det(/-C2) > w^.w^ ^^^^'' ^"^^^^ I I n M^) ^"^*^^) I ^.^ Mze) 

~ Im(a) lm{a) Im(a) Im(a) Im(a) Im(a) 

^ + (B.8) 

im(aj im(a) 

We now find an upper bound on the perturbation term e(z). To this end, we have 
< W2 < •^+dmax/^^ and 0<V2< •^+dLax/^z- Recalling dH]), we also have 

h _ {^21 < 1 + ^""^^ <^max Qmax I ^ I 



Using the same arguments as in the proof of Proposition I6.2l -(i) (involving this time the 
tightness of the measures associated with the Stieltjes transforms ^Tr DR and ^TyDR) 
yields: 

Im(a) - ' l^^^^^l - n'^^l^l^^^"- ^ ' Ini(a) " ^'d^in ' 

for every z G C+ U C^, where 77, are positive constants, and $ and ^, nice polynomials. 
Finally, we can state that there exist nice polynomials $ and vj/ such that: 

By continuity of det(/ — C2{z)) at any point of the open real negative axis, this inequality 
is true for any z G C — R+. Denote by £n the set: 

£n = l^zeC-R+, (^1^ < 1/2 

If z G £m then det(/ — C2) is readily lower-bounded by the quantity stated in Proposition 
K2\(u). 

By considering inequalities (jB.6P and ()B.7P and by possibly modifying the polynomials $ 
and we have 1 — M2 > and 1 — U2 > for z G £"„. The proof of proposition 16.21 - (ii) is 
completed. 
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